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Abstract — The purpose of designing a controller for a 
teleoperation system is achieving stability and high efficiency 
in presence of factors such as time delay, noise, system 
disturbance and modeling errors. In this article seven new 
schemes for teleoperation systems is suggested to reduce the 
error of tracking between the master and the slave systems. In 
the first and second structures, optimal control is applied to the 
master and in the third structure it is applied to the slave. In 
the forth structure, optimal controller has been designed in 
both the master and slave subsystems and by a suitable 
combination of the output signals of both controllers and 
exerting it to slave, it has been tried to create the best 
performance with regard to tracking. In the 5th structure, 
optimal controller is applied to both the master and slave 
systems and the output of each controller is then applied to it's 
own system. In the structures 6 and 7, using optimal control 
and wave transformation, it has been tried to reduce the 
tracking error between the master and the slave systems which 
is considered as the cost function in this article. In these two 
structures, optimal controller guarantees the efficiency of the 
system and wave variables method has been used to guarantee 
the stability of the system for constant time delay. In this 
method, the appropriate selection of weighting factors of cost 
function and wave impedance plays a decisive role in the 
improvement of system performance and reduction of tracking 
error between the master and the slave. 
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I. Introduction 

A teleoperation system enables human operator to 
implement given task in a remote manner or enhance his/her 
capability to handle both the micro and the micro worlds. A 
typical teleoperation system consists of a local master 
manipulator (master site) and a remotely located slave 
manipulator (slave site). The human operator controls the 
local master manipulator to drive the slave one to implement 
a given task remotely. The system must be completely 
"transparent" so that the human operator can feel as if she/he 
were able to manipulate the remote environment directly. In 
general, in the design of teleoperation systems there is a trade 
off between high transparency and sufficient stability 



margins [1]. The main control strategies proposed for 
bilateral teleoperation systems with constant delay include 
the delay compensation proposed by Anderson and Spong 
[2], which assures the asymptotic stability by Niemeyer and 
Slotine [3], where wave transformations are used to keep the 
passivity of the communication channel. Remote compliance 
control proposed by Kim et al. [4], where the remote 
compliance smoothes the interaction between the slave and 
the environment, but the transparency is poor for delays 
greater than 1 s. The internet has time- varying characteristics 
(for example, the amplitude of time delay), therefore it is 
useful to use information-measured in-line of communication 
channel to design the control structure, as in [5,6]. Where 
these delay compensations are based on the wave variables 
modifying any main parameter (apparent damping or gain, 
respectively) depending on the time delay. The control 
structures based on wave variables keep the passivity with 
delay impedance, but at the expense of reducing system 
transparency. The paper is organized as follows: in section 2 
wave variables method is introduced and in section 3 the 
optimal control design used in this article is presented. In 
section 4 simulations of different structures for different 
weighting factors is demonstrated and the validity of this 
scheme is stablished. In section 5 the results are discussed. 

II. WAVE VARIABLES METHOD 

Wave variables method presents a modification or 
extension to the theory of passivity, which creates robustness 
to arbitrary time delay. This method is also closely related to 
the scattering and small gain theories [7]. 

In this sense, we define U to denote the forward or right 
moving wave, while V denotes the backward or left moving 
wave. 

The wave variables (u,v) can be computed from the 
standard power variables (x,f) by the following 
transformation. 



v = b -^,v- 



42b 



bX-F 
42b 



(1) 



Where b is the characteristic wave impedance and may be a 
positive constant or a symmetric positive definite matrix. 
Fig. 1 shows the basic wave transformation. 
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Fig. 1 : the basic wave transformation 

using the velocity X and left moving wave V as the 
inputs, the transformation will determine the force F and 
right moving wave U. All other combinations are also 
possible. 



F = bX-yl2bV , U = -V + V2bX (2) 

Fig. 2 shows wave variables transformation used in this 
article. 
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Fig. 2: wave based communication, by transforming velocity-force variables 
(subscripts s and m denote slave and master respectively) 

In this configuration, both velocity and force variables 
are first transformed into wave variables before transmission. 
Upon arrival information (velocity and force) is extracted 
from the slave variables, in accordance with the convention 
used by Niemeyer [8], right moving waves are defined as 
positive and left moving waves are negative. The full details 
about wave variables method is described in [7]. 

III. OPTIMAL CONTROL DESIGN IN TRACKING 
PROBLEMS 

Continues time linear system is presented in the state 
space like [9]: 



X(t) = AX (t)+Bu(t) 
Cost function that should be minimized: 

t f 
J = ^({X(t)-r(t)fQ{X(t)-r(t)) 

+ -U \t))Ru{t\dt 



(3) 



(4) 



The final t M fixed. H and Q are real symmetric positive 

semi-definite matrices and R is real symmetric positive 
definite matrix. It is supposed that the controls and positions 
are not bounded and x(t f ) is free. The acceptable control 

u(t) that causes the equation system of (3) to follow the 
desired route of r(t) and minimize the cost function is as 
follows: 



u(t) = -R~ 1 B T k(t)X(t)-R~ 1 B T s{t)] 

= F(t)r(t)+V(t) 



(5) 



In which F(t) is the feedback reinforcement matrix and 
V(t) is the command signal. It should be noted that V(t) 
depends on the system parameters and the desired route(r(t)). 
In fact V(t) depends on the future amounts of the reference 
signal. 

K (t) and S (t) in the above equations are: 

k(t) = -K(t)A - A T K(t) -Q+ K(t)BR V x(t) (6) 

S(t) = -[A T -K(t)BR~ 1 B T ]s(t)+Qr(t) (7) 

And boundary conditions are: 

K(t f ) = H (8) 

4)="4) w 

Since in this problem the value of H is zero, as a result 
the boundary conditions in the equations (8) and (9) are zero. 

To achieve u(t), we should integrate equation (6) and (7) 
from t to t f and using boundary conditions we will 

find u(t) . In the first and second structure we design an 

optimal control in the Master system and in the third 
structure we design an optimal control in the slave system. 
The reference signal for the optimal controller in the first 
structure is X s (t-r) and for the second structure is X m (t) 
and the desired signal for the optimal controller in the third 
structure is x m (t - r). If we place the reference signals of the 

controllers in equation (6) and (7) we can find u(t) from 

equation (5). Figure 3 (4 th structure) shows a new structure of 
bilateral teleoperation system with time delay in the 
transmission channel using combinational optimal controller. 
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Fig. 3: structure of teleoperation system using combinational optimal 
controller (4th structure). 



In this structure optimal controller has been designed in 
both the master and slave subsystems and the controlling 
signals of each controller have been multiplied by a certain 
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coefficient and are added together and the resulted 
controlling signal is applied to the slave in a way that it can 
follow the master in the least time. In this scheme we have 
considered only one degree of freedom for the slave and the 
master in the teleoperated system. The dynamic model of an 
element with a dof is as below: 



J6(t) + b0(t) + mglsm 0{t) = u(t) 



(10) 



J is the inertia of the element, M is the mass of robot 
manipulator, g is the gravity acceleration, 1 is the length of 
robot manipulator, #(j)is the rotate angle, b is the viscous 

friction coefficient, and u(t) is the control torque applied. The 
simplified linear model is: 



J0(t) + bd(t) = u(t) 



(11) 



Considering the position (Xj(/)= o(t)) and the velocity 
(x 2 (t)= 6(t)) as state variables, the state space representation 
of the master and slave is: 



Xm 2 (t) 

x Sl (t)' 

Xs 2 (t) 



'Xm x (t)' 
Xm 2 (t) 

x Sl (t)' 

Xs 2 (t) 



u m (t) 



u s (t) 



y m 



(0 = [i o] 



: (0=[i o 



Xm^(t) 

Xm 2 (t) 

ix Sl (t) 

[Xs 2 (t) 



(12) 



(13) 



(14) 



(15) 



The structure of the other matrices that appear in the 
model, in this case, is the next: 



Kn 



-[k, 



in k m2 \ 



Rm ~ Vm\ r mi\ 



(16) 



(17) 



In which K m is the feedback matrix of the master state. 
R m defines the interaction between the slave-master. It 
allows modeling the force reflection to the master. 



In the teleoperation systems, it is 
account the remote environment 
manipulator performs a task can 
environment. This interaction must 
design of the generated reaction forces 
simplification of the Kelvin model of 
by means of the stiffness ( k ) and a 

The reaction force that acts against the 

f s {t)=k e e s {t) + b e e s (t) 



necessary to take into 
. When the slave 
interacting with the 
be considered in the 
. We have employed a 
the environment only 
viscous friction (b e ). 
slave is given by: 

(18) 



In this way the reaction force, f s (t), will oppose to the 

slave control signal, modeling that the slave is interacting 
with the environment. To consider force feedback from the 
slave to the master, the structure of the matrix R m must be 
[10]: 



K(t)=[r ml r m2 ]=[k f k e k f b e ] 



Where k * is the force feedback. 



(19) 



The master control signal, that is u m (t) , as shown in 
figure (3) is: 



u m (t) = K m (t)x m (t)+ R m (t)x s (t -T)+ F m (t) 



(20) 



If we place the above control signal in the state space 
representation of the master we have: 



x m (t)= (A m + B m K m )x m (t) 
+ B m R m x s (t-T)+B m F m (t) 



(21) 



And for the slave system, the control signal should be 
calculated so that the slave system can optimize the cost function 
of equation (4) by tracking the desired route of r(t). the control 
signal of the slave system is: 



u s (t)=F(t)x s (t)+V(t) 



(22) 



If we place the above control signal in the state space 
representation of the slave we have: 

*,(*)= A s*s(t)+ B s [F(t)x s (t)+v(t)] (23) 

In which v(t) depends on the slave system parameters 
and the reference signal (x m (t - r)) 

F(t) and v(t) in the above equation are calculated from 
Riccati equations in equation (6 ) and (7). 

In the 5th structure, the same as 4th structure, optimal 
controller is applied to both the master and slave systems but 
the output of each controller is then applied to it's own 
system. The reference signal for the optimal controller in the 
slave subsystem is X m (t-r) and the desired signal for the 
optimal controller in the master system is X s (t-r). If we 
place the reference signals of the master and the slave 
systems in equation (6) and (7) we can find u{t) . 

The following figure shows the 6th designed structure for 
teleoperation systems using wave variables and optimal 
control methods. 
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Fig. 4: Structure of teleoperation system using combinational optimal 
controller and wave transfirmation 

In this structure, the position and velocity variables on 
the master side are transformed into wave variables. Then 
they pass the transmission channel and are transformed to 
position and velocity variables on the slave side again. The 
returning signal on the slave side are also first transformed 
into wave variables, then passing the transmission channel, 
wave variables are transformed into position and velocity 
variables on the master side. The optimal controller is used 
on both master and slave sides. The optimal control on the 
master side applies the controlling signal to the master so 
that the master output follows the delayed slave output. The 
optimal control on the slave side functions so that the slave 
output signal follows the delayed master output. The output 
of each controller is applied to its own system so that the 
tracking error between the master and slave is minimized. 

In the 7th structure designed for teleoperation systems 
using optimal control methods and wave variables, an 
optimal controller is used on the slave side in a way that the 
slave output follows the delayed master output. 

IV. SIMULATION RESULTS 

The considered parameters of the master and the slave 
are respectively: 

J m =1.5kgm 2 ,b m =ll Nm f m m =28A25kg, l m =0Am 

rad/ 

/ s 



Therefore regarding the given parameters, state matrices 
of the slave and the master are as follows: 



A m = 




1 

L5 



-11 



[1 0] 



[0 


1 1 





-15 

2 


r°~ 




1 


,C S = 


_2_ 





[1 0] 

Given the equation (19), the value of the matrix R m is: 
R m =[k f k e k f b e ] = [l0 O.l] 

Considering that the poles of the master are placed in the 
position -15 of the s plane, the following control parameters 
are obtained: 

KAK* U=[-249 -33 

In this paper all of the structures are simulated for 
different weighting factors. The best weighting factors are 
achieved through trial and error. 

The simulation results are demonstrated in figures 5 to 
10. 
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Fig. 5: .simulation of master- slave behaviors (a: input; b: position; c: 
position error; d: time delay)for the first structure (T=1.5s) 



J s =2kgm 2 , b s =15 



Nm 
rad / 



m =1.5kg, l s =2m 



To consider the interaction between the slave and remote 
environment, we have considered a stiffness ( k e ) and a 
viscous friction (b ) as below: 



b^lJ^ . k e =100Nm/ ad 



rad/ 



Considering that the force feedback gain is 0. 1 
* /= 0.1 
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Fig. 6: simulation of master- slave behaviors (a: input; b: position; c: 

position error; d: time delay) for the second structure with variable time 

delay 



© 2010 ACEEE 
DOI:01.ijepe.01.01.09 



49 



4cACEEE 



ACEEE International Journal on Electrical and Power Engineering, Vol. 1, No. 1, Jan 2010 




Time(sec) 



Time (sec) 



Fig. 7: simulation of master- slave behaviors (a: input; b: position; c: 
position error; d: time delay) for the third structure 
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Fig. 8: simulation of master- slave behaviors (a: input; b: position; c: 

position error; d: time delay) for the5 th structure using combinational 

optimal controller in master and slave syb systems 



Fig. 10: simulation of master- slave behaviors (a: input; b: 
position; c: position error; d: time delay) for the 7 th structure using optimal 
controller and wave variable transformation 



V. Conclusion 

In this paper, we have proposed the design of a bilateral 
controller for teleoperation systems with variable time delay. 
These controllers ensure the passivity, transparency and high 
efficiency of the system and take into account the 
communication delays and the interaction between the slave 
and the environment. In these structures, the optimal 
controller receives delayed position signals of master and 
slave and compares them to the reference signals and applies 
control signals to the master and slave in a way that the 
tracking error between them is minimized and they can track 
each other in the least possible time. For stabilizing the 
system for variable time delay, wave variable method has 
been used. The simulation results have shown the superiority 
of the augmented technique to the commonly used wave- 
variable method. Simulation results have shown excellent 
performance in motion tracking and good response in force 
tracking. The model uses a state space representation of the 
teleoperation systems considering all the possible interaction 
that can appear in the operator-master-slave-environment set. 
Comparative experiments demonstrate the validity of the 
proposed controls method and their excellent performance of 
motion/force tracking. 
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Fig. 9: simulation of master- slave behaviors (a: input; b: 
position; c: position error; d: time delay) for the 6 th structure using optimal 
controller and wave variable transformation 
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